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ABSTRACT 


yy 


Part  I  showed  that  (a)  the  type  of  anomaly  was  not  critical  in 
computing  deflection  of  the  vertical  if  the  Molodenski  theory  was  used 
and  complete  gravity  knowledge  existed,  (b)  what  was  critical  was 
extrapolation  of  gravity  values  between  observation  points,  (c)  the 
complete  Bouguer  anomaly  with  geologic  corrections  provided  the  anomaly 
which  allowed  the  most  reliable  extrapolations,  and  (d)  the  normal 
complete  Bouguer  anomaly  provided  the  most  readily  derived  values  for 
computation  of  deflections.  In  this  part  of  the  final  report.  Part  II, 
the  procedure  outlined  in  Part  I  is  applied  to  the  problem  of  inter¬ 
polation  of  deflections  of  the  vertical. 

In  the  Test  Phase  &f  the  eon tract  three  astro-geodetic  deflection 
stations  in  the  Rocky  Mountain  area  of  the  western  United  States^  wespe 
Aj-hiH  These  were  deflection  stations  102,  I05,  and  II6  from 

USC&GS  Special  Publication  No.  229*  Using  the  adopted  procedure,  the 

i, 

%  deflection  component  was  interpolated  between  stations  105  and  II6 
to  obtain  an  interpolated^  value  at  station  102.  A  comparison  of 
this  interpolated  \  value  with  the  astro-geodetic^'  value  showed  that 
the  two  differed  by  only  ,22^,  Since  the  accuracy  of  the  astro- 
geodetic  values  probably  do  not  exceed  +.2^,  the  results  were  considered 


highly  satisfactory 

in^the  Application  Phase 


» t vv  s iv e 


deflection  interpolation 


woo  oftrrtod  Qut  in  an  area  in  the  A3ps  chosen  by  AGIO.  In  this  case 
three  deflection  stations  were  chosen  from  a  deflection  of  the  verti¬ 
cal  map  provided  by  AGIO.  These  were  designated  stations  1,  2,  and  3 
for  convenience.  Using  the  adopted  procedure,  the  X‘l®^'l®ctlon  com¬ 
ponent  was  Interpolated  between  stations  1  and  3  to  obtain  a  value  at 


station  2.  As  there  were  no  individual  deflection  values  given,  the 


f 


astro-geodetic  deflections  were  estimated  from  a  deflection  component 

/  ■  ‘ 

contour  map  with  a  2.50’'  contour  interval.  Because  of  the  manner  of 
obtaining  the  astro-geodetic  deflections,  they  are  considered  to  have 
an  accuracy  of  the  order  of  +  .5"* 

/ 

i'OL’  this  reason  the  difference  between  the  computed  and  observed 
deflection  components,  .59”  >  has  an  uncertainty  of  +1.0"^  Since 
the  difference  again  lies  within  the  range  of  error  inherent  in  the 
data  the  results  are  considered  entirely  satisfactory. 
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SECTION  1 


GENERAL  COMMENTS 

In  Part  1  of  the  final  report  on  work  completed  under  contract 
AF23( 601) -4009  the  final  formula  for  deflection  of  the  vertical  compu¬ 
tation  was  given  as 

(0  =  -1.  r  (ig'  +  2fllcah)  ~  iti  +  ABp  (1-1) 

where  cp  is  the  gravimetric  deflection  in  the  angular  direction  9 

o 

R  =  mean  radius  of  the  earth 
Yq  =  theoretical  gravity  at  deflection  station 
Ag*  =  complete  Bouguer  anomaly 

k  =  universal  gravitational  constant 

a  =  density  used  in  computing  Bouguer  anomalies 
h  =  elevation  of  surface 
S  (i)|  )  =  Stokes  function 

^  -  central  angle  between  deflection  station  and  incremental 
surface  element  dw 

AQp  =  correction  factor  dependent  on  elevation  near  deflection 
station  dw  =  sin  ijf  diJfdA 

A  =  azimuth  angle  from  computation  point 
The  formula  for  A  6p  can  be  stated 

3  '5 

ASp  =  L-H  J  .  I  cos  I  cos  (r^,l)  drdA  (1-2) 

r  r 

where  r  ==  distance  from  deflection  station  to  incremental  surface 
element,  cos (r^,  l)  =  angle  between  the  direction  in  which  the  deflec¬ 
tion  is  being  computed  and  direction  of  the  surface  element,  dw. 

E  =  h  -  Sq  =  elevation  difference  between  station  and  surface  element 
The  aim  in  the  Test  and  Application  Phases  of  the  contract  was 


to  test  formulae  (l-l)  and  (l-2)  in  making  interpolations. of . deflections 
of  the  vertical.  In  making  the  test  deflection  of  the  vertical  inter¬ 
polations,  optimum  use  was  made  of  both  astro-geodetic  and  gravimetric 
information  for  computing  a  deflection  at  a  point  lying  between  sta¬ 
tions  where  astro-geodetic  deflections  were  available.  If  one  had 
conrplete  knowledge  of  the  gravity  field  and  carried  out  the  integrals 
of  (l-l)  and  (1-2)  around  the  world  to  determine  a  gravimetric  deflec¬ 
tion  at  an  astro-geodetic  station,  the  gravimetric  deflection  would 
be  identical  to  the  astro-geodetic  deflection  —  assuming  the  astro- 
geodetic  deflection  was  acciarate  and  was  referred  to  the  same  reference 
ellipsoid  as  that  used  in  computing  the  gravity  anomalies#  The  major 
part  of  the  deflection  components  at  a  station  is  due  to  the  values  of 
the  integrated  quantities  of  equations  (l-l)  and  (l-2)  within  a  short 
distance  of  the  deflection  station.  This  contribution  of  the  near  area 
varies  irregularly  from  point  to  point#  The  effect  of  the  gravity 
field  at  greater  distances,  on  the  other  hand,  has  less  effect  on  the 
deflection  and  is  a  smoothly  varying  function# 

If  we  label  the  astro-geodetic  deflection  as  cp  and  the  gravl- 

a 

metric  deflection  as  cp  we  can  write 

g 


'fa  =  =  'fg„  'fg, 

n 


(1-3) 


where  cp  represents  the  part  of  the  gravimetric  deflection  obtained 
Sn 

by  integration  of  (l-l)  and  (l-2)  over  the  area  near  the 
deflection  station,  and  lying  within  a  circle  having  a  radius 
r^^  described  about  the  station, 

and  cp  represents  the  part  of  the  gravimetric  deflection  obtained  by 

Of 

integration  of  (l-l)  and  (1-2)  over  the  area  lying  at  a  dis¬ 
tance  greater  than  r^  from  the  deflection  station.  The  value 

of  r^^  will  depend  upon  how  much  area  it  is  desirable  to  retain  in  the 
"near"  part  of  the  deflection,  cpg^#  In  general,  the  value  of  r^  will 
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be  controlled  by  the  smoothness  of  the  gravity  field  and  the  distance 
over  which  deflection  interpolation  is  to  be  carried  out. 

To  obtain  a  deflection  at  on  intermediate  station  lying  between 
two  astro-geodetic  stations,  cp  is  computed  at  each  astro-geodetic 
station  by  numerically  evaluating  the  integrals  of  (l-l)  and  (1-2) 
between  zero  and  r^. 

Then  at  each  astro-geodetic  deflection  station  we  can  compute 


tPf 


<Pa  -  <pg_ 


n-it) 


Since  is  a  smoothly  varying  function,  one  can  compute  the  value  of 
at  the  intermediate  station  by  linear  interpolation  between  the  two 
astro- geodetic  stations.  Then  ?p  is  computed  at  the  intermediate  sta- 
tion  in  the  same  way  as  it  was  at  the  astro -geodetic  stations,  added 
to  the  interpolated  cp^  value,  and  a  complete  deflection  of  the  vertical 
which  we- shall  designate,  obtained. 

In  many  cases  the  orientation  and  origin  of  the  ellipsoid  to  v/hich 
the  astro-geodetic  deflections  are  referenced  differs  from  the  Inter¬ 
national  Ellipsoid  to  which  the  gravity  values  are  referenced.  So  long 
as  one  desires  the  intermediate  deflection  to  be  referenced  to  the  same 
ellipsoid  as  the  astro- geodetic  deflections,  this  presents  no  problem 
in  carrying  out  the  interpolation.  In  the  case  of  a  difference  in  ref¬ 
erence  ellipsoid,  equation  (l-3)  becomes 


‘(’a  “  '<’g  “  ‘Pg  'fgj 


(1-5) 


where  Acp^  represents  the  angle  between  the  normals  to  the  two  ellip¬ 
soids  at  the  point.  Then  at  each  astro-geodetic  station  we  have 
tpf  =  'I’a  -  ‘•’g  “  '(’f  + 

Since  Acp  is  a  smoothly  varying  function,  we  can  interpolate  to  get 
e 

at  the  intermediate  station  and  then  compute  cpg  there  to  again 
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arrive  at  a  cp  value  at  the  intermediate  station  which  is  referenced 
to  the  astro- geodetic  ellipsoid. 

Rice  has  developed  a  method  of  numerically  computing  the  value 
of  an  integral  of  the  form 


(1-7) 


where  F  can  he  any  type  of  gravity  anomaly.  This  method  can  be  used 
to  evaluate  the  integral  of  equation  (l-l). 

In  Rice's  method  a  template  such  as  indicated  schematically  below 
is  used. 


The  template  is  so  set  up  that  for  each  of  the  compartments 
beyond  the  central  circle  one  need  only  multiply  the  average  gravity 
anomaly  in  mgls  of  that  compartment  by  the  cosine  of  the  angle  between 
the  direction  of  the  center  of  that  compartment  and  the  direction  in 
which  the  deflection  is  being  computed  and  then  by  .001  to  give  the 
deflection  contribution  of  that  compartment  in  seconds. 

The  effect  of  the  central  circle  is  then  computed  by  using  values 
of  gravity  on  the  circle's  perimeter  to  estimate  the  gravity  gradient 
within  the  circle.  The  details  of  the  template  and  central  circle 
computation  as  used  by  Rice  are  given  in  Heiskanen  and  Vening  Heines z 
(1958)  and  are  well  knovm  to  most  geodesists.  They  v^ill  not  be  re¬ 
peated  here. 
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In  the  present  problem  the  Rice  template  is  used  twice  to  compute 
the  effect  of  the  area  outside  the  central  circle. 

First  the  integral 

-R_  ig'  ^  sin  ijldildA  (1-3) 

5«Yq  Jo 

is  evaluated  using  Rice*s  template  with  average  values  of  Ag'  being 
p1  -Tor  eaoh  compartment. 

The  integral 


fo  4^  2,tkch  M  sin  $d(,dA 


(1-9) 


I,  360 
-  R_  P  ^ 

"  o 

is  then  evaluated  using  Rice’s  template  in  the  following  manner.  A 
normal  template  calculation  is  carried  out  using  the  average  elevation 
for  each  compartment  rather  than  a  gravity  anomaly.  The  final  result 
is  then  multiplied  by  2itka  to  conveit  from  feet  to  mgls  and  thus  give 
the  correct  deflection  contribution  of  the  integral. 

A  central  circle  calculation  is  then  carried  out  to  determine 
the  contribution  of  the  central  circle  to  each  of  the  integrals  of 
(1-8)  and  (1-9). 


As  shOT-m  in  Part  I  of  this  final  report,  for  computation  of  the 

quantity  A  0  ,  one  may  sum  the  effect  of  a  number  of  annular  compart- 

P 

ments, either  the  compartments  of  the  Rice  template  or  others.  The 
effect  of  a  single  compartment  is  given  by 


)  =  ~  AA  cos(r  .  _i,)-  3h^  (-K  -  -K)} 


(1- 


where  A  A  is  the  angular  width  of  the  compartment 

and  r^j^  and  r^  are  the  inner  and  outer  radii  of  the  compartme 
One  could,  of  course,  establish  any  values  one  desired  for  aA 


10) 


and  for  the  various  radii  increments.  If,  for  example,  we  use 
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AA  =  15°  «=  .2617  radians,  then  equation  (l-lO)  becomes 

A(A0  )  =*  +.  694  7  X  10"^  cos  -  ”’2) 

P  r2 

-2,42  X  10- ®  cos(r^,l)h^(-i5  .  -^)  (1-11) 

*^1  "^2 

vliei-e  h  is  given  in  units  of  thousands  of  feet  and  r;j_  and  Pg  are  given 
in  km.  If  h  is  in  km  the  constants  in  equation  (l-ll)  become 

+  .2117x10*^  and  -  .8l86xl0"^ 

_3 

Ue  shall  call  the  term  containing  h  the  secondary  correction  and 
the  term  containing  h^  the  tertiary  correction.  Let  us  see  how  these 
quantities  vary  with  r.  We  can  irrlte  the  secondaiy  correction  as 

+  .6947  X  10“^  cos(r^,  *  F(rj^,r2)  cos(r^,l)h^ 

where  F^(r^,r2)  »  ,6947  x  10' ^  - -^)  (1-12) 

r ,  ro 

Then  the  following  table  gives  the  value  of  F(r2^,  r^)  for  various 
values  of  r^  and  r2  in  kms  where  h  must  be  given  in  units  of  thousands 
of  feet. 


“l 

”^2 

^2 

) 

km 

km 

.2 

.5 

14.5337 

X 

10"  ^ 

.5 

.75 

1.5436 

X 

10-2 

.  75 

1.1 

.6614 

CM 

1 

O 

1.  1 

1.55 

.  2341 

X 

10-  2 

1.55 

2.  IS 

.  1452 

X 

'10-2 

2.  IG 

3.07 

.0709 

X 

10-2 

3.07 

4.32 

.0520 

X 

10-  2 

The  tertiary  correction  can  be  written 

-2,4202  X  10'®  cos  (rg,l)h^(--j - i^)  =  G(rj,r2)  cos(r^,l)h^ 

ri 

where  G(rj^,r2)  =•  -2,4202  x  10-®(~j,  -  -i,,) 

ri  r2 


(1-  13) 
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Then  the  values  of  G(rjLi  r^)  for'-lrarious  r2  values  are 
given  belov  where  again  it  is  assumed  and  r2  are  in  kms  and  h  in 
thoiAsands  of  feet. 


'1 

'2 

G(ri,  tj) 

km 

km 

.2 

.5 

.1474  X  10'^ 

.5 

.75 

.0031  X  10'^ 

.75 

1.1 

.0006  X  10"^ 

.1 

1.55 

.00012  X  10"^ 

The  above  tables  show  that  the  secondary  correction  can  be  com-  ! 

puted  sufficiently  accurately  by  integrating  only  over  an  area  within 

5  km  of  the  deflection  station  in  the  normal  case.  The  tertiary  correc-  i 

tion  normally  need  only  be  integrated  over  an  area  within  1  km  of  the  I 

I 

deflection  station.  Indeed,  unless  the  elevation  changes  involved  are  i 

I 

rather  large,  the  A  0  p  correction  can  be  ignored.  I 

I 

The  solution  given  in  equation  (l-lO),  as  was  the  case  for  the  | 

Stokes*  equation,  cannot  be  used  within  a  small  central  circle  around  j 

the  deflection  station.  No  attempt  was  made  to  obtain  a  central  ring 

solution  for  equation  (l-lO).  It  was  felt  that  the  assumption  that  the  i 

t 

area  within  .2  km  (about  600  ft.)  of  a  deflection  station  is  sufficiently 

level  to  introduce  no  significant  error  was  reasonable.  | 

I 

The  first  question  which  naturally  arises  in  carrying  out  the 
deflection  interpolation  is  the  question  of  how  far  it  is  necessary  to 
carry  out  the  numerical  integrations  of  equations  (1-8)  and  (1-9)  so 

I  " 

that  the  or  value  obtained  can  be  linearly  interpolated  between 
the  astro- geodetic  stations.  The  answer  to  this  question  depends  of 
course  upon  the  distance  between  astro-geodetic  stations.  Molodensky 
et  al  (1980)  suggest  that  as  a  minimum  the  integration  be  carried  out 
to  a  distance  from  each  station  equal  to  the  separation  of  the  astro- 
geodetic  stations.  Of  course  theoretically  the  further  the  integration 
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is  carried  out  the  more  nearly  linear  will  .l>e  the  interpolated 
quantity.  In  both  the  Test  and  Application  Phases  of  this  study 
the  integration  was  carried  out  to  a  distance  approximately  equal  to 
the  separation  of  the  astro -geode tic  deflection  stations.  This 
appeared  to  be  adequate  to  obtain  an  answer  whose  accuracy  was  con¬ 
sistent  with  the  accuracy  of  the  gravity  and  deflection  data  used  in 
the  computations  and  comparisons.  Where  modem  astro-geodetic  deflec¬ 
tions  are  available  with  an  accuracy  of  .1"  or  better  and  a  reasonably 
dense  net  of  high  quality  gravity  stations  are  available,  it  might  be 
desirable  to  extend  the  integration  over  a  circular  area  whose  radius 
was  about  twice  the  distance  between  deflection  stations  in  order  to 
aim  at  an  interpolation  accuracy  of  between  .1"  and  .2". 

In  evaluating  the  integrals  of  equations  (1-8)  and  (1-9)  using  the 
Rice  template  and  inner  circle. method,  it  is  always  necessary  to  decide 
the  radius  of  the  inner  circle.  Such  a  decision  must  be  made  in  each 
case  on  the  basis  of  the  smoothness  of  the  variation  of  the  quantity 
being  integrated  near  the  deflection  station. 
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SECTION  2 

TEST  PHASE  -  METHOD  OF  PROCEDURE 

In  order  to  test  the  accuracy  of  the  selected  deflection  interpolation 
procedure,  application  of  the  theory  to  actual  data  was  undertaken.  Three 
astro- geodetic  deflection  stations  were  selected  in  the  western  United 
States  to  use  in  carrying  out  the  test.  These  were  deflection  stations  102, 
105,  and  Il6  from  USC&GS  Special  Publication  No.  229  (l9^l). 

The  test  was  carried  out  as  follows:  the  meridional  deflection,  §  , 
was  considered  to  be  known  at  stations  116  and  105  and  the  interpolation 
procedure  performed  to  obtain  a  computed  deflection  at  station  102,  This 
computed  deflection  was  then  compared  with  the  astro -geode tic  deflection 
at  station  102  as  a  measure  of  the  accuracy  of  the  deflection  interpola¬ 
tion  procedure. 

Steps  followed  in  deri\’lng  the  interpolated  deflection  values,  were 
as'  follows: 

1.  The  area  within  65*9  of  each  deflection  station  (through 
Rice  Ring  3T)  was  chosen  as  the  near  area  over  which  equations  (l-l) 
and  (l~2)  were  to  be  numerically  integrated  to  obtain  g  values, 

•y 

2.  For  the  area  between  4.32  km  and  65*9  km,  elevation  contour 
maps,  scale  1:250,000,  were  used  to  obtain  average  elevations  for  the 
compartments  of  a  Rice  template.  These  average  elevations  were  used 
to  evaluate  the  integral  of  equation  (1-9)  in  this  region. 

3.  To  obtain  the  summation  of  (1-9)  over  the  area  within  4.32  km 
of  the  stations,  the  following  procedure  was  used.  For  deflection 
stations  ll6  and  102,  an  elevation  map  whose  scale  was  1:62,500  was 
used  with  a  Rice  template  to  pick  average  elevations  for  evaluating  the 
integral  between  1,099  km  and  4.32  km  from  the  deflection  station. 

The  Rice  inner  circle  procedure  employing  the  three  gradient  method 
was  used  to  compute  the  effect  of  the  area  within  1.099  km  of  these 
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two  stations.  In  the  case  of  deflection  station  I05,  an  inner  circle 
procedure  employing  the  single  gradient  method  was  used  for  the  entire 
area  within  4.32  km  of  the  station  since  the  topography  was  irelatively 
nmcx)th  near  this  station. 

4.  Using  the  elevation  map  of  scale  1; 62, average  elevations 

were  determined  for  use  in  computing  A  9  at  station  II6  employing 

P 

equation  (l-lO).  The  details  of  the  size  of  sectors  over  which  averag¬ 
ing  was  carried  out  are  discussed  with  the  results.  Preliminary  calcu¬ 
lations  showed  that  the  A  9  correction  would  be  negligible  at  stations 

P 

105  and  102,  thus  it  was  not  computed  for  these  stations. 

5.  Simple  Bouguer  anomalies  were  plotted  on  a  transparent  paper 
overlay  for  the  1:250,000  elevation  contour  maps.  Terrain  corrections 
were  then  computed  using  the  procedure  developed  under  contract  AF  23 
(601) -3739  and  complete  Bouguer  anomalies  obtained. 

6.  Using  geophysical  and  geologic  information  to  control  the  inter¬ 
polation  between  observation  points,  a  Complete  Bouguer  anomaly  contour 
map  was  prepared  on  the  scale  1:250,000. 

7.  Using  a  Rice  template,  average  anomalies  for  compartments  were 
determined  for  the  area  between  4.32  km  and  65.9  km  and  used  to  evaluate 
the  integral  of  equation  (I-8)  in  this  area.  The  Rice  inner  ring  pro¬ 
cedure  was  used  to  evaluate  (I-8)  over  the  area  \dthin  4.32  km  of  each 
station. 

8.  The  final  computation  of  the  part  of  the  meridianal  deflection, 
component, §  due  to  the  gravity  field  of  the  near  area  (within  65.9 
km  of  a  station)  was  made  for  each  station  using  the  results  obtained 
in  the  preceding  steps. 

9.  At  stations  13j6  and  105,  the  deflection  contribution  of  the  near 
area  was  subtracted  from  the  astro-geodetic  deflections.  The  result. 


was  the  deflection  contribution  of  the  gravity  field  of  the  area 
beyond  65.9  tan  plus  the  effect  of  the  differences  between  the  gravi¬ 
metric  and  astro-geodetic  ellipsoids. 

10.  Using  a  simple  linear  interpolation  between  Stations  II6  and 

105,  the  value  was  determined  at  Station  102.  This  interpolated 

value  was  then  added  to  the  previously  computed  deflection  effect  of 

the  near  area  §  at  Station  102  to  give  a  computed  meridonal  deflection 
Sn 

component  at  Station  102.  This  was  then  compared  with  the  astro-geo¬ 
detic  value,  §a,  to  test  the  accuracy  of  the  interpolation  procedure. 

The  above  steps  give  in  outline  form  the  manner  in  which  the 
deflection  interpolation  procedure  was  carried  out  for  the  Test  Phase 
of  the  Contract.  The  paragraphs  below  give  additional  details  concern¬ 
ing  the  computations  indicated  in  the  above  outline. 

Elevation  Contour  Map  Selection 

The  first  step  in  evaluating  the  integral  of  equation  (1-9)  vas 
the  selection  of  elevation  contour  maps  from  which  average  elevations 
could  be  determined  for  the  template  compartments.  The  choice  of  map 
scale  was  controlled  by  the  amount  and  quality  of  the  gravity  data 
present.  An  elevation  contour  map  should  be  chosen  which  would  allow 
the  average  elevations  of  the  compartments  of  the  template  to  be  deter¬ 
mined  with  an  accuracy  such  that  the  quantities  2nkoh  computed  from  them 
would  be  about  equal  in  accuracy  to  the  average  Complete  Bouguer  anomalies 
determined  for  the  same  compartments.  It  was  felt  that  the  compartment 
averages  of  Complete  Bouguer  anomaly  could  be  determined  to  an  accuracy 
of  about  +  3  mgls.  Thus  an  elevation  accuracy  of  about  +  100  ft.  was 
desired.  For  the  area  outside  Rice  ring  21  (outer  radius  4.32  tan),  the 
1:250,000  series  of  Transverse  Mercator  Projection  maps  prepared  by  the 
Army  Map  Service  and  sold  by  the  U.  S.  Geological  Survey  was  considered 
adequate.  For  the  area  within  4.32  km  of  a  deflection  station,  a  larger 
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scale  map  was  desired.  For  the  present  problem  the  1:62,500  special 
tcpogrcphic  quadrjingle  nt-pti  of  the  U.  S,  Ecological  Survey  '.rere  chosen. 

Transverse  Mercator  Projection  maps  were  used  throughout  the  com¬ 
putations.  Transverse  Mercator  Projection  maps  cannot,  of  course,  be 
used  in  carrying  out  world -wide  computations,  but  they  are  accurate 
enou^  to  be  used  for  interpolation  computations  where  only  areas  of 
1  to  2  degrees  square  need  be  used  in  the  integration. 

The  manner  in  which  the  average  elevations  of  the  compartments 
were  determined  varied  depending  upon  the  size  of  the  compartment  and 
the  ruggedness  of  the  topography.  For  the  smallest  compartments  and 
in  level  areas  the  average  compartment  elevation  vras  estimated  directly 
by  visual  inspection.  For  the  larger  compartments,  particularly  where 
the  topography  was  rugged,  the  compartment  was  divided  into  a  number 
of  subsections  —  up  to  nine  subsections  in  the  most  extreme  case. 

The  average  elevations  of  these  subsections  were  then  determined  and 
themselves  averaged  to  obtain  the  average  elevation  for  a  compartment. 
Preparation  of  Bouguer  Anomaly  Map 

To  obtain  the  Complete  Bouguer  anomaly  contour  map  necessary  to 
use  for  determination  of  compartment  averages  utilizing  a  Rice  template, 
the  following  procedure  was  chosen.  The  simple  Bouguer  anomalies  were 
first  plotted  on  a  transparent  tracing  paper  overlay  on  the  1:250,000 
Transverse  Mercator  elevation  contour  maps .  Then,  using  a  procedure 
developed  under  contract  AF  23(601) -3789>  terrain  corrections  were  made 
at  all  stations  where  it  was  felt  that  the  correction  exceed  1.5  mgls. 

The  terrain  correction  procedure  is  not  described  here  as  it  will  be 
found  in  the  final  report  to  contract  AF  23(601) -3789»  Once  the  com¬ 
plete  Bouguer  anomalies  had  been  computed  and  plotted,  the  contouring 
of  the  complete  Boug’jer  anomaly  map  was  undertaken.  It  was  at  this  point 
that  geologic  and  geophysical  knowledge  were  interjected.  The  geologic 
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and  geophysical  knowledge  was  used  to  control  the  interpolation  (l.e., 
conto\n*ing)  between  points  of  observation.  The  primary  utilization  of 
geologic  knowledge  was  in  contouring  the  gravity  data  in  the  Rocky 
Moiintain  Front  Range  of  Colorado  and  in  the  area  of  transition  from 
the  Denver- Julesburg  Basin  to  the  Front  Range  Uplift.  To  gain  insight 
into  the  form  the  gravity  field  should  have  in  passing  from  the  basin 
to  the  uplift,  a  geologic  density  section  based  on  available  geologic 
knowledge  was  prepared  and  a  gravitational  profile  computed.  Utiliz¬ 
ing  these  results,  the  available  gravity  data,  and  knowledge  of  the 
location  of  the  fault  boundary  between  the  basin  and  the  uplift,  the 
contouring  in  this  area  was  carried  out.  On  the  western  side  of  the 
Front  Range  uplift  not  enough  information  was  found  to  justify  numer¬ 
ical  computation  of  the  geologic  effect.  The  best  that  could  be  done 
was  to  utilize  the  available  geologic  knowledge  of  structural  rela¬ 
tions,  approximate  thickness  of  sediments,  and  location  of  recent 
Cenozoic  intrusives  to  control  the  contouring.  The  geologic  control 
utilized  in  the  contouring  process  although  largely  non-quantitatlve 
rather  than  quantitative  proved  highly  successful  as  witnessed  by 


the  results  obtained 
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SECTION  3 

TEST  PHASE  -  RESULTS 


In  Tables  1  through  3  of  the  Appendix  are  presented  the  average 
Cmplete  Bouguer  anomalies  and  average  elevations  for  each  compartment 
of  Rice  circles  22  through  37  for  the  three  deflection  stations  used  in 
the  Test  Phase.  In  each  case  the  results  for  similar  compartments  of 
each  ring  were  summed  and  each  sum  multiplied  by  the  appropriate  cosine. 
The  sum  times  cosine  results  were  in  turn  summed  to  give  a  single 
result  for  each  of  the  two  quantities  at  each  station.  The  final  Complete 
Bouguer  anomaly  result  at  each  station  was  multiplied  by  .001  to  give 
the  Complete  Bouguer  anomaly  component  in  seconds  of  arc.  The  final 
average  elevation  result  was  multiplied  first  by  .03406  to  convert  to 
mgls  and  then  by  .001  to  obtain  the  Bouguet  correction  contribution  to 
the  deflection  component.  The  results  are  listed  below. 

Station  ll6 

Bouguer  Correction  Component  +  1.290" 

Complete  Bouguer  Anomaly  Component  +  .829" 

Station  IO5 

Bouguer  Correction  Component  -  3*659" 

Complete  Bouguer  Anomaly  Component  -  .084" 

Station  102 

Bouguer  Correction  Component  -  3*450" 

Complete  Bouguer  Anomaly  Component  +  2.012" 

For  the  circular  area  between  the  deflection  stations  and  Rice 
ring  22  (O  to  4.32  km)  the  effects  of  the  two  components  were  computed 
using  such  combinations  of  template  simmation  and  central  circle  compu¬ 
tation  as  were  considered  necessary  to  obtain  adequate  accuracy. 

The  Complete  Bouguer  Anomaly  Component  was  computed  at  each  sta¬ 
tion  by  a  Rice  central  circle  computation.  The  results  of  these  com¬ 
putations  are  given  belovr. 


At  Stations  ll6  and  105  the  Complete  Bouguer  anomaly  contour  lines 


run  almost  directly  north-south.  Thus  at  these  stations  the  Complete 

Bouguer  anomaly  component  for  the  central  circle  vels  found  to  he  zero. 

For  Station  102,  the  single  gradient  method  (See  Heiskanen  and 

Vening  Heines z,  1958)  was  used  to  get  the  Complete  Bouguer  anomaly 

component  of  the  inner  ring.  Using  the  formula 

(A  g  -  Ag  ) 

Ag  =  .  105  r  — § - 


and  the  values  r  =*  4 . 3  2  km 
o 


Ax 


ve  get 


Ax  »  8.64  km  Ag  =  -  219  mgls 
=  .105  =  0,42" 


=  -  227  m 


The  elevation  contour  map  which  had  to  be  integrated  to  obtain  the 
Bouguer  correction  component  had  a  more  complex  contour  pattern  than  the 
Complete  Bouguer  anomaly  component  and  required  the  use  of  a  combination 
of  template  and  central  circle  computation  to  evaluate  the  effect  of  the 
eirea  \7ithin  4.32  km  of  each  statj.on.  For  stations  ll6  and  102,  1:62,500 
scale  maps  were  used  with  a  Rice  Template  to  obtain  the  Bouguer  correc¬ 
tion  component  for  Rice  zones  l4  through  21.  The  results  are  given  in 
Tables  4  and  5  of  the  Appendix.  The  deflection  contributions  resulting 
from  this  template  summation  were 

Station  ll6  +  1.50" 

Station  102  -  .32" 


The  small  inner  circle  of  1.1  km  radius  was  evaluated  at  each  of  these 
stations  using  a  three  gradient  method.  The  results  given  below  were 


computed  using  the  formula 

-  ^gn> 


A§"  =  0.105r 


=  0.003576  r 


Ax 


-  En) 

Ax 


southward  elevation 
northward  elevation 


A?"  =  deflection  contribution  of  central  circle 
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Station  ll6 

Tq  =  1.1  km  12,100 


FOR  CENTRAL  INTERVAL 

Eg  -  =  13,150  -  12,100  =  1050 

LX  =  2.2  )sm 

and  .  =  (0.105)  (.03406)  =  (0.105)(.034o6)(525) 

Ll*^  =  ( 55. 12)(. 03406)  =  +  1.88" 

FOR  RIGHT  INTERVAL 

Eg  -  =  (12,400  -  12,600)  =  -  200 

LX  =  1.5556 

=  (0.105)(l.l)  =  -  .506" 

1.5556 

FOR  LEFT  INTERVAL 


Eg  -  En  =  (12,900  -  13,250)  =  -  350 


AX  =  1.5556 

=  (0.105)  (1*1)(»03406)(-  350)  =  .  ,505" 

^  1.5556 

1 

ig  AVERAGE  =  ^  (Al^  +  =  1  (l,88  -  ) 

AS  AVERAGE  =  i  (I.88  -  .JO)  =  |  {I.I8)  =  .59" 


Station  102 

=  1.1  km 


5,950 
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FOR  CENTRAL  INTERVAL 

E  -  E  =  (5950  -  6020)  =  -  70 
s  u 

AX  =  2.2  km 

A§"c  =  (0.105)(1.1)(. 03406) 

A§"c  =  -  .125 

FOR  RIGHT  INTERVAL 

E^  -  E„  =  0 
s  n 

A§"j,  .  0 
FOR  LEFT  INTERVAL 

Eg  -  Ejj  =  (5970  -  6000)-  =  '-  30 

AX  =  1.5556 

A§"l  =  (0.105)(1.1)  '  =  -  -OTS 

Ag  AVERAGE  =  5  ( A§o  =  =  |  (.  .125  - 

Al  AVERAGE  =  |  (-  .125  -  -038)  =  5  (-  .163)  =  -  .08I 

For  Station  IO5  the  Bouguer  correction  component  vas  computed  using 
the  single  gradient  method  for  the  area  within  4.32  km  of  the  deflection 
station  since  the  topogr?,phy  varied  in  a  linear  manner  near  this  station. 
The  results  are  given  helo^r. 


6700 


6400 

r  =  4.32  km  •  AX  =  8.64  km 

ThenAg"  =  (  .105)(  .03406)  '  ^7°°) 

Ag"  =  -  (.io5)(.034o6)(150) 

Ag"  =  -  (.003576)(150) 

Ag"  =  -  .536" 
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As  pointed  out  in  Seciion  1,  very  large  changes  in  elevation  are 
required  very  near  a  station  in  order  for  the  correction  factorAO  ^  to 
he  sufficiently  large  to  require  computation.  For  this  reason  only  at 
Station  ll6  did  the  component  have  a  significant  contribution  to  the 
final  deflection.  For  this  station  A 8p  uas  computed  using  equation  (l-2) 
and  a  template  with  the  radii. -of  .the  rings  being  as  follows: 


Inner 

radius 

Outer  radius 

Ring  A 

.2 

km 

.5 

km 

Ring  B 

.5 

km 

.75 

km 

Ring  C 

.75  km 

1.1 

km 

Ring  D 

1.1 

km 

2.0 

km 

Ring  E 

2.0 

km 

3.0 

km 

Ring  F 

3.0 

km 

h.O 

km 

Rings  A,  B,  and  C  were  divided  into  annular  sections  of  10 
degrees  as  in  the  Rice  template.  Rings  D,  E,  and  F  were  divided  into 
annular  sections  of  15  degrees.  The  area  -v^ithin  .2  km  was  assumed  to 
have  zero  contribution. 

As  stated  in  Part  I  of  the  final  report,  average  elevations  obtained 
using  a  Rice  template  would  normally  be  used  so  that  average  elevation 
determinations  already  available  from  determination  of  the  Bouguer  cor¬ 
rection  component  could  be  utilized.  In  the  present  case,  the  computa¬ 
tion  of  the  A^  component  at  Station  ll6  had  already  been  accomplished 
using  the  template  divisions  described  above  before  it  was  decided  that 
the  results  of  the  Rice  template  computations  would  have  been  satisfac¬ 
tory.  Thus,  this  solution  was  retained  since  recomputing  would  not  have 
contributed  to  any  increase  in  accuracy.  The  results  of  the  A9^  compu¬ 
tations  are  given  in  Tables  6  and  J  in  the  Appendix,  The  contribution 
of  this  term  to  the  deflection  at  Station  ll6  was  -  ,665". 
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V/ith  the  results  given  above  ve  are  nov  in  position  to  carry  out 
the  deflection  Interpolation.  Below  is  given  a  summary  of  the  above 
results  at  each  station. 


Station  ll6 


Bouguer  Correction  Component  (4.32  km  to  65.9  hm) 
Bouguer  Correction  Component  (I.099  km  to  4.32  km) 
Bouguer  Correction  Component  (O  to  1.099  km) 

TOTAL  Bouguer  Correction  Component 
Complete  Bouguer  Anomaly  Component  (4.32  km  to  65.9 
Complete  Bouguer  Anomaly  Component  (O  to  4.32  km) 
TOTAL  Complete  Bouguer  Anomaly  Component 
AQp  Component 

TOTAL  value  of  f 


+  1.290" 
+  1.500" 
+  .590" 
+  3.380" 
km)  +  .829" 
.000" 
+  .829" 
-  .814" 

+  3.395" 


Station  105 


Bouguer  Correction  Component  (4.32  km  to  65.9  km)  “  3.659" 

Bouguer  Correction  Component  (O  to  4.32  km)  -  .536" 

TOTAL  Bouguer  Correction  Component  -  4.195" 

Complete  Bouguer  Anomaly  Component  (4.32  km  to  65.9  km)  -  .084" 
Complete  Bouguer  Anomaly  Component  (O  to  4.32  km)  .000') 

TOTAL  Complete  Bouguer  Anomaly  Component  -  .084'.' 


A  9p  Component 
TOTAL  vlaue  of 
Station  102 


«  .000" 
-  4.279 


Bouguer  Correction  Component  (4.32  km  to  65.9  km)  -  3.450" 

Bouguer  Correction  Component  (1.099  km  to' 4. 32  km)  -  .320" 

Bouguer  Correction  Component  (O  to  1.099  km)  -  ,08l" 

TOTAL  Bouguer  Correction  Component  -  3.85I" 
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Complete  B’ouguer  Anomaly  Component  (4.32  tan  to  65.9  km)  +  2.012*' 
Complete  Bouguer  Anomaly  Component  (0  to  4.32  km)  +  .420" 

TOTAL  Complete  Bouguer  Anomaly  Component  -«•  2.432" 


Component 
TOTAL  value  of  § 


.000" 


^n 


-  1.419" 


The  computation  of  §  _  +  A§  gives  the  following  results  for  the 

Sf  6 

two  stations,  II6  and  IO5,  at  which  astro-geodetic  deflections  were 
assvimed  known: 

Station  II6 


Astro-geodetic  deflection 
Effect  of  near  region 

Sn 

-  Sg,  =  Sg,  - 

Station  IO5 

Astro-geodetic  deflection  S 


Effect  of  near  region  § 
a  Sn  gf  a 


Sn 


+  1.35" 

+  3.40" 

-  2.05" 

-  5.34" 

-  4.28" 

-  1.06 


Figure  3-1  shows  the  linear  interpolation  carried  out  to  obtain 

a  value  of  §_  +  A§  for  Station  102.  The  value  obtained  for  Station 

of  e 

102  is 

=  -  ^-69 

But  the  value  of  computed  for  Station  102  was  -  1.42".  Thus 
on 

the  interpolated  deflection  for  Station  102  is 
=  -  1.69"  -  1.42"  =  -  3.11" 

The  astro-geodetic  deflection  for  Station  102  is  -  2.89".  Thus  the 
interpolated  deflection  value  and  the  observed  deflection  value  differ 
only  by  .22". 
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SECTION  k 

APPLICATION  PHASE  -  METHOD  OF  PROCEDURE 

In  order  to  test  the  applicability  of  the  deflection  interpolation 
procedure  developed  under  this  contract  to  general  usage  in  elII  areas 
of  the  world  the  procedure  was  used  to  interpolate  deflections  of  the 
vertical  in  the  European  Alps,  Gravity  data  and  astro-geodetic  deflec¬ 
tions  for  this  phase  of  the  study  were  provided  by  ACIC,  The  gravity  . 
data  consisted  of  2000  gravity  stations  located  over  eleven  1°  X  1° 
squares  in  the  European  Alps.  The  astro -geode tic  deflection  stations 
where  selected  in  conjunction  with  ACIC  from  stations  shovm  on  a  con¬ 
tour  map  of  T]  deflection  components.  These  astro- geodetic  stations  were 
for  convenience  labeled  #1,  7^2,  and  /r3»  Their  location  as  determined  by 
scaling  from  the  contour  maps  were  as  follows: 

station#!  Lat,  4?°  22.5'H  Lon^.  8°  to.5'E 

Station  #2  Lat.  47°  22.5’W  Long.  9°’01.5'E 

Station  #3  Lat.  47°  22.5'W  Long.  9°  39.0'E 

Because  of  the  large  number  of  gravity  stations  and  the  excellent 
elevation  data  available  it  was  hoped  that  the  application  of  the  present 
deflection  interpolation  procedure  in  the  Alpine  area  would  provide  an 
excellent  test  of  the  method.  However,  the  test  did  not  prove  as  con¬ 
clusive  as  was  first  hoped.  There  were  two  primary  reasons  for  this. 
First,  the  form  in  which  the  astro- geodetic  information  was  available 
placed  a  limit  on  the  extent  to  which  the  accuracy  could  be  tested.  The 
values  of  the  astro -geode tic  deflection  components  could  be  determined 
from  the  contour  maps  only  to  an  accuracy  of  about  +  ,5".  Thus  it  was 
only  possible  to  determine  if  the  accuracy  of  the  interpolated  deflection 
was  within  the  range  of  the  cumulative  uncertainty  which  was  approximately 
+  1.0".  Also,  because  of  the  scale  of  the  deflection  contour  map. 
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1  inch  =  70  hm,  uncertainty  of  about  +  1.0’  exists  in  the  location 
of  the  stations  as  scaled  from  these  maps.  This  introduces  uncertain¬ 
ties  into  the  computed  components  which  are  estimated  to  be  in  the 
range  of  .1"  to  .3”. 

A  second  problem  encountered  was  in  the  distriubtion  and  accuracy 
of  the  gravity  data.  It  was  originally  anticipated  that  the  2000 
stations  would  lie  approximately  within  a  3°  X  3^  square  or  nine  1*^  X  1^ 
squares.  Instead  the  data  was  located  within  eleven  1°  X  1°  squares. 
Moreover,  nearly  half  the  stations,  about  900,  were  located  within  one 
15*  X  15'  square  area.  After  plotting  of  the  data  and  computation  of 
the  terrain  corrections  an  attempt  was  made  to  contour  the  data.  At 
this  time,  t\70  more  problems  were  encountered.  First,  it  was  found  that 
at  many  stations  there  were  two  gravity  anomaly  values  given  which 
differed  by  2.0  mgls.  A  check  of  latitude,  longitude,  and  elevation 
data  showed  that  this  discrepancy  evidently  resulted  from  the  listing 
of  the  same  survey  twice  with  different  base  values  used  for  each 
listing.  One  set  of  values  was  chosen  and  the  other  set  rejected  on  the 
basis  of  comparisons  with  other  stations. 

Second,  it  was  found  that  in  some  places  closely  spaced  stations 
differed  by  several  tens  of  milligals.  Also  in  certain  areas  where  only 
a  few  stations  occurred,  it  was  found  that  contouring  using  the  available 
data  produced  a  gravity  anomaly  contour  map  which  had  no  apparent  rela¬ 
tion  to  geology  and,  in  fact,  appeared  to  crosscut  geologic  structure. 
Such  a  gravity  pattern  appeared  to  be  highly  unlikely.  The  cause  of 
the  above  discrepencies  was  traced  to  certain  series  of  very  old  pre-1900 
pendulum  data.  This  data  was  in  error  by  amounts  up  to  50  mgls.  This 
series  of  stations  wets  discarded  and  the  contour  map  redraim.  However, 
because  of  the  factors  listed  above,  there  were  left  available  for  con¬ 
touring  the  complete  Bouguer  anomaly  map  even  fewer  useful  data  points 
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than  had  been  available  for  the  Test  Phase. 

The  above  factors  although  disappointing  in  one  sense  do  tend  to 
emphasize  certain  useful  facts.  Foremost  is  the  fact  that  vhere  obser¬ 
vational  data  is  suspect,  the  employment  of  geologic  knowledge  can  lead 
to  criteria  for  acceptance  or  rejection  of  the  data.  Second  is  the  fact 
that  even  a  very  small  amount  of  accurate  data  when  judiciously  handled 
will  give  reasonably  good  deflection  results. 

Since  the  details  of  procedure  for  the  Application  Phase  of  the 
contract  were  similar  to  those  of  the  Test  Phase  they  vrLll  not  be 
brought  out  here.  Certain  special  details  do,  however,  merit  mention. 

The  area  used  in  computation  of  the  T[„  components  is  somewhat 

Sn 

smaller  than  would  normally  be  used  if  precise  answers  were  desired. 

That  is,  if  one  were  attempting  ,1"  accuracy  idth  plentiful  gravity  data, 
integration  should  be  carried  out  over  several  additional  Rice  rings  than 
were  used  in  the  present  case.  Considering  the  accuracy  of  the  astro- 
geodetic  data,  however,  such  an  extension  in  the  integration  area  would 
not  be  expected  to  yield  any  additional  useful  information  in  the  present 
case.  In  fact,  as  a  check,  the  integration  was  extended  through  an 
additional  Rice  ring  (ring  37)  after  all  the  tables  included  in  this 
report  had  been  completed. 

These  extensions  changed  the  actual  values  of  by  amoiints  in 
the  range  .1"  to  .2"  and  altered  the  interpolated  value  by  less  than 
.1".  This  indicated  that  the  range  of  change  to  be  expected  by  exten¬ 
sion  of  the  area  of  integration  was  so  small  as  to  be  of  q^uestionable 
meaning  considering  the  sparse  gravity  control.  In  any  case,  the  uncer¬ 
tainty  in  the  astro- geodetic  deflections  precludes  any  meaningful  inter¬ 
pretation  of  the  effect  of  including  the  additional  area  on  the  accuracy 
of  interpolated  deflections. 
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The  elevation  contour  maps  used  for  the  determination  of  average 
compartment  elevation  in  areas  beyond  4.32  km  from  deflection  stations 
were  Army  Map  Service  1:250,000  Transverse  Mercator  Projection  maps, 
series  1959*  For  the  area  idLthin  4.32  km  of  the  deflection  stations, 
German  1:50,000  scale  elevation  contour  maps  were  used  to  estimate 
average  elevations. 


SECTION  5 

APPLICATION  PHASE  -  RESULTS 

The  results  of  the  Application  Phase  of  the  Contract  are 
outlined  here  with  the  various  tables  given  in  the  Appendix. 

The  Bouguer  Correction  Component  for  all  three  stations  was  com¬ 
puted  for  the  area  between  4.321mi  and  55 *661011  using  the  1:250,000 
scale  maps.  The  Complete  Bouguer  Anomaly  Component  was  computed 
using  the  Complete  Bouguer  anomaly  contour  map  prepared  on  the  1:250, 
000  scale.  The  average  compartment  elevation  values  used  to  compute 
the  Bouguer  correction  component  are  given  in  Table  0,  9^  and  10  in 
the  Appendix.  The  average  Complete  Bouguer  anomalies  for  the  com¬ 
partments  are  given  in  Tables  11,  12,  and  13  of  the  Appendix. 

The  final  results  for  the  above  summations  bet^/een  4.321an  and 
55.66km  at  the  three  deflection  stations  are  given  below: 

Station  #1 


Bouguer  Correction  Component  -2.3^’ 

Complete  Bouguer  Anomaly  Component  +2.47" 

Station  #2 

Bouguer  Correction  Component  -1.62" 

Complete  Bouguer  Anomaly  Component  +1.07" 

Station  #3 

Bouguer  Correction  Component  5.05" 

Complete  Bouguer  Anomaly  Component  +3.11" 


For  station  #1  a  three  gradient  computation  was  carried  out  for  the 
Bouguer  Correction  Component  and  a  single  gradient  computation  for 
the  Complete  Bouguer  Anomaly  Component  within  the  4.32km  radius  central 


circle 
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The  results  were: 

Complete  Bouguer  Anomaly  Component 

-89  'o  =  ‘^•32  %  =  -  ^ 

_ /  ATI'  =  .0525  C(-a6)  -  (.89)]  =  +  .16" 


Bouguer  Correction  Component 


r  =  4.32 

o 

central  interval  ATf'  =  -  ,59" 
north  interval  ATI"  =  -  .67" 
south  interval  AT]"  s  +  .17” 
Average  ATf'  =  -  .42" 


For  station  #2  a  single  gradient  central  circle  computation  was 
carried  out  to  obtain  the  Complete  Bouguer  Anomaly  Component  of  the 
inner  circle.  The  Bouguer  Correction  Component  was  obtained  by  using 
a  template  summation  between  1.304lna  and  4.321aii  and  a  three  gradient 
central  circle  calculation  for  the  area  within  l,099toa  of  the  sta¬ 
tion.  The  average  elevation  results  of  the  template  calculation  are 
given  in  Table  14  of  the  Appendix.  The  Bouguer  Correction  Component 
of  the  deflection  for  the  area  between  1.304km  computed  from  these 
average  elevations  was  +  .606".  The  results  of  the  two  central  circle 
computations  were: 

Complete  Bouguer  Anomaly  Component 
ro  =  4.321nn  Agv^  =-101  Age  =-102 
ATf’  =  .0525  il-lOl  -  (-102)] 

A'H"  =  .052" 

-101 1  }  -102 
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» 


Bouguer  Correction  Conponent 
Tq  =  1.304km 


north  interval  ATI" 

south  interval  A'l]” 

Average  AT]" 


.461" 


For  station  #3  a  single  gradient  central  circle  computation  was 
carried  out  to  obtain  the  effect  of  the  4.32km  inner  circle  for  the 
Complete  Bouguer  Anomaly  Component.  The  Bouguer  Correction  Compo¬ 
nent  was  foimd  to  he  zero.  The  results  for  the  Complete  Bouguer 
Anomaly  Component  are  given  below. 

Complete  Bouguer  Anomaly  Component 


Tq  =  4.321aii  gw  =  -  120;  gg  =  -  12? 
ill"  =  .0525  [(-120)  -  (-127)3  t 


ATI"  =  .357" 


-127 


The  correction  component  A9  p  was  found  to  be  negligable  at  all 
of  the  stations. 

To  carry  out  the  interpolation  the  quantity  was  ccmputed  at 


the  three  stations.  Summarizing  the  results  we  have  below; 

Station  1 

Bouguer  Correction  Component  (4.32km  to  55.66km)  -2.34" 

Bouguer  Correction  Component  (O  to  4.32km)  -  .42" 

TOTAL  Bouguer  Correction  Component  -2.76" 

Complete  Bouguer  Anomaly  Component  (4.32km  to  55.66km)  +2.47" 
Complete  Bouguer  Anomaly  Component  (0  to  4.32km)  +  .16" 

TOTAL  Complete  Bouguer  Correction  Congjonent  +2.63" 

A6p  Component  0.00" 

TOTAL  Tlgp^  for  Station  1  -  .13" 
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Station  2 

Bouguer  Correction  Component  (4.32k]n  to  55 *661011)  -1.62" 

Bouguer  Correction  Component  (l.304lnn  to  4.321nii)  +  .82" 

Bouguer  Correction  Component  (O.  to  1.304km)  +  .46" 

TOTAL  Bouguer  Correction  Component  -  *34" 


Complete  Bouguer  Anomaly  Component  (4.321an  to  55  *661011)  +1,87" 
Complete  Bouguer  Anomaly  Component  (l.099ton  to  4.321aii)  +  .05" 


TOTAL  Complete  Bouguer  Anomaly  Component  +1.92" 

A6p  Component  O.OO" 

TOTAL  for  Station  2  +  1.53'’ 

Station  3 

Bouguer  Correction  Con^ionent  (4.321aii  to  55*66lan)  -5*85" 

Bouguer  Correction  Component  (O  to  4.321a3i)  0*00" 

TOTAL  Bouguer  Correction  Component  -5*85 

Complete  Bouguer  Anomaly  Component  (4.321an  to  55*66lan)  +3*11 
Complete  Bouguer  Anomaly  Component  (O  to  4.321nn)  *36 

TOTAL  Bouguer  Correction  Component  +3*4? 

A 0p  Component  0.00 

TOTAL  value  for  Station  3  -2.38 


Because  of  the  fact  mentioned  previously  that  the  T1  values  were 
picked  from  a  contour  map  the  values  determined  for  the  astro-geodetic 
deflections  have  considerable  imcertainty.  The  values  estimated  for 
the  astro-geodetic  deflections  were; 

Station  1  Ti  =  -  I.50"  ±  .5" 

a 

station  2  Tla  =  -  2,00"  ±  .5" 

Station  3  Tla  =  -  ±  *5" 

Using  the  above  values  and  the  Tgjj  values  computed  previously 
the  following  VEilues  for  Tlj  were  determined  at  stations  1  and  3» 
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station  1  Tlf  =  Ha  -  1160  =  '  ^•3'^"  -  *5 

.  Station's  HJ  =  Ha  -  H6n  =  -  5.62"  +  .5 

As  shown  in  Figure  5*'l  interpolation  of  the  above  values  leads 
to  an  interpolated  value  at  Station  2  of  =  -  2 •91"  ±  .5*  this 
is  added  to  the  value  of  Tlgj^  computed  at  Station  2  we  obtain  an  inter¬ 
polated  value  of  T]  =  Tl^  +  Tlg^^  =  -  1,1*1  ±  .5*  The  difference  between 
the  interpolated  and  the  astro -geodetic  deflection  is  therefore 
T|a-T|=-,59  +1.0".  The  relation  between  the  interpolated  and  the 
astro-geodetic  deflection  component  is  illustrated  in  a  slightly  dif¬ 
ferent  manner  in  Figure  5-1  where  the  computed  value  of  T]  ^  at  Station  2 
is  compared  with  the  interpolated  value. 
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SECTION  6 

SUMMARY  AND  CONCLUSIONS 

In  both  the  Test  and  Application  Phases  of  this  study,  'the  re¬ 
sults  obtained  with  the  proposed  deflection  interpolation  method  were 
excellent.  In  each  case,  the  agreement  between  interpolated  and  observed 
deflection  components  was  close  enough  that  it  was  within  the  limits  of 
sensitivity  of  the  test  data.  The  sensitivity  of  the  Test  Phase  was 
controlled  by  the  accuracy  of  the  astro- geodetic  def lection^ components . 
The  probable  error  claimed  for  the  astro-geodetic  deflections  was  +  .1" 
to  +  .2".  The  results  of  the  deflection  computations  of,  for  example. 
Rice  (1952),  hov/ever,  show  that  the  accuracy  of  astro- geodetic  deflec¬ 
tions  is  often  much  less.  Since  the  astro-geodetic  deflections  used 
were  the  results  of  measurements  made  prior  to  I920  a  probable  error 
of  ±  .2"  would  more  than  likely  be  optimistic.  Under  these  conditions 
agreement  of  the  interpolated  and  the  astro-geodetic  deflections  to 
vrLthin  ,22"  can  certainly  be  considered  as  within  the  limits  of  possi¬ 
ble  error  in  the  astro-geodetic  values. 

As  mentioned  previously  the  accuracy  of  the  astro-geodetic  data 
used  in  the  Application  Phase  was  such  that  only  errors  in  the  inter¬ 
polated  deflections  greater  than  1.0"  could  be  clearly  determined.  The 
astro- geodetic  minus  interpolated  T]  deflection  component  difference  of 
.59"  ±  1.0"  is  V7ell  within  the  limits  of  sensitivity  of  the  data. 

The  accuracy  of  the  results  obtained  with  any  gravimetric  deflec¬ 
tion  interpolation  procedure  must  be  dependent  upon  the  amount  and 
quality  of  the  gravity  data  available.  In  this  respect  the  data 
available  for  both  the  Test  and  Application  Phases  leaves  much  to  be 
desired.  Indeed,  the  results  would  appear  to  be  better  than  one  should 
normally  expect  from  the  gravity  coverage  available.  Unfortunately 
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extensive  gravity  data  is  not  likely  to  be  available  in  most  rugged 
mountainous  areas  in  the  near  future.  Especially  gravity  data  is  not 
likely  to  be  available  in  close  spacing  near  astro-geodetic  deflection 
stations  unless  special,  expensive  efforts  are  made  to  obtain  it.  It 
1b y  therefore,  important  to  know  if  the  apparent  success  obtained  vith 
the  chosen  deflection  interpolation  method  using  a  minimum  of  gravity 
data  is  the  result  of  a  basic  characteristic  of  the  method  or  singly 
a  result  of  luck. 

The  answer  would  seem  to  be  a  little  of  both.  Although  the  writers 
would  be  the  first  to  admit  that  the  degree  of  accuracy  obtained  is 
unlikely  to  be  always  possible  using  comparable  amounts  of  gravity 
data  the  deflection  interpolation  procedure  chosen  here  does  appear 
to  allow  optimum  utilization  of  available  data  and  to  render  unneces¬ 
sary  the  carrying  out  of  extremely  close  spaced  s\.irveys  near  deflection 
stations. 

The  i*eason  satisfactory  results  can  be  obtained  •v/ith  the  chosen 
method  -vdthout  an  excess  of  gravity  data  can  be  seen  from  the  nature 
of  the  equations  and  has  been  mentioned  in  the  discussing  of  Part  I 
of  the  final  report.  The  gravitational  effect  of  the  topography  is 
separated  from  the  rest  of  the  gravitational  field  and  interpolated 
between  gravity  stations  by  using  a  topographic  map.  In  few  if  any 
instances  would  it  be  possible  to  obtain  sufficient  gravity  data  to 
make  direct  interpolation  of  the  gravitational  effect  of  topography 
from  the  gravity  maps  more  accurate  than  that  obtained  using  elevation 
contour  maps  and  an  assumed  density  based  on  geologic  knov/ledge.  In 
the  present  problems  the  standard  density  of  2.67gni/cc  has  been  used 
in  the  computations  as  being  as  good  as  any  other.  In  areas  such  as 
the  Gulf  Coast  where  geologic  knowledge  indicates  that  the  density  of 
the  majority  of  the  topographic  features  is  not  2.67gm/cc^the  correct 
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density  should  be  used. 

With  the  topography  accounted  for  one  is  interested  only  in  ob¬ 
taining  sufficient  data  to  accurately  interpolate  Complete  Bouguer 
anomalies  between  observations.  The  Complete  Bouguer  anomaly  is  always 
a  smoothly  varying  function  and  thus  requires  considerably  fewer  sta¬ 
tions  for  its  definition.  By  using  geologic  and  geophysical  informa¬ 
tion  much  of  the  variability  in  the  Complete  Bouguer  anomaly  can  be 
understood  and  therefore  better  estimates  of  the  variation  between 
observation  points  made.  Although  how  many  or  how  few  observation 
points  are  needed  will  ultimately  depend  upon  the  local  situation,  it 
is  hard  to  imagine  that  there  are  many  cases  when  a  10km  grid  of  sta¬ 
tions  with  perhaps  5  additional  stations  \7ithin  10km  of  the  deflection 
station  would  not  prove  completely  adequate.  Of  course,  one  is  un¬ 
likely  to  often  have  gravity  on  a  10km  grid  but  the  above  estimation 
is  meant  simply  to  show  the  order  of  station  density  at  which  one  should 
aim.  The  point  is  that  the  addition  of  station  density  around  a  de¬ 
flection  station  so  that  there  are  more  than  say  20  stations  at  a  maxi¬ 
mum  TTithin  10km  of  a  deflection  station  is  unlikely  to  lead  to  a  very 
significant  increase  in  the  accuracy  with  which  the  gravity  field  can 
be  derived.  Attempts  to  obtain  gravity  stations  at  1km  spacing  appear 
to  be  completely  superfilous. 

The  results  obtained  were  very  encouraging  and  indicate  the  method 
used  is  a  good  one.  To  establish  the  method’s  actual  degroc^cf  reliability, 
a  more  extensive  test  should  be  run  involving  the  use  of  sets  of  three 
stations  not  in  a  line  to  interpolate  both  deflection  components  to  a 
fourth  station  lying  within  the  triangle  of  which  the  three  stations  are 
vertices.  Eventually  this  type  of  interpolation  might  be  used  not  only 
to  establish  intermediate  delfection  stations  but  to  test  the  accuracy 
of  astro- geodetic  deflection  stations  and  thus  improve  the  astro-geodetic 
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def lection  network  by  discovery  and  rejection  of  Incorrect  values. 

As  the  primary  purpose  of  the  solutions  obtained  under  the 
present  contract  was  to  test  the  accuracy  of  the  deflection  Interpo¬ 
lation  procedure  using  Bouguer  anomalies,  It  Is  essential  that  a  stand¬ 
ard  computational  method  such  as  the  Rice  template  summation  be  used  In 
obtaining  the  deflection  component  results. 

Recommendations  on  Procedure 

During  the  course  of  this  work  serious  consideration  was  given  to 
developing  a  computer  technique  for  carrying  out  the  computations  along 
lines  of  the  method  utilized  by  Nagy  (I962),  However,  a  rather  large 
programming  job  with  a  considerable  amount  of  trial  and  error  in  order 
to  arrive  at  an  optimum  procedure  is  involved.  In  addition  elevation 
and  Bouguer  anomaly  data  would  need  to  be  averaged  on  an  entirely 
different  basis. 

As  useful  results  could  not  be  obtained  during  the  contract 
duration  and  as  this  was  not  an  integral  part  of  the  contract  the 
attempt  was  reluctantly  abandoned.  It  is  believed,  however,  that 
except  for  the  area  within  about  I5  km  of  the  deflection  station  such 
a  computer  technique  will  ultimately  be  required  for  large  scale 
problems.  The  computer  method  will  become  peuftlcularly  Important 
in  the  near  future  when  very  accurate  5*  x  5*  average  elevations 
become  available  for  large  areas  such  as  the  continental  United 
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EXPIAIIATIOK  OF  TABLES 


Final  results  can  be  obtained  for  Bouguer  Correction  Component 
and  Complete  Bouguer  Anomaly  Component  from  the  tables  in  this 
Appendix  as  follows. 

Test  Phase 

Bouguer  Correction  Component  =  -  ( .03406)  ( .001)  SIM  X  COS 
Complete  Bouguer  Anomaly  Component  =  +  (.001)  SUII  X  COS 

Application  Phase 

Bouguer  Correction  Component  ■  -  ( .111?) ( .001)  SUli  X  SIN 
Complete  Bouguer  Anomaly  Component  =  +  ( .  OOl)  SUI'I  X  SIN 

The  change  in  sign  between  the  Bouguer  Correction  Component  and 
the  Complete  Bouguer  Anomaly  Component  results  from  the  fact  that  in 
the  Tables  the  Complete  Bouguer  Anomalies  have  been  listed  as  posi¬ 
tive  numbers  to  avoid  needless  repetition  of  the  minus  sign  tlirough- 
out  the  tables  although  the  Complete  Bouguer  Anomalies  used  are  nega¬ 
tive.  This  is  taken  into  account  with  a  sign  change  in  the  summation 
equations  above. 
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Table  4 

DEFLECTION  STATION  102 


BOUGUER  CORRECTION  COMPONENT 


14 

RICE 

•RISGS-. 

14-21 

15 

16 

17 

18 

19 

20 

21 

SUM 

( 

SUM  X  COsj 

6025 

6030 

6060 

6060 

6060 

6100 

6110 

6130 

48575 

48390 

6012 

6025 

6040 

6050 

6060 

6100 

6110 

6130 

48527 

46872 

6025 

6050 

6050 

6075 

6090 

6100 

6100 

6180 

48670 

44109 

6033 

6060 

6075 

6100 

6060 

6150 

6160 

6210 

48848 

40016 

6050 

6060 

6080 

6100 

6150 

6200 

6190 

6200 

49030 

34669 

6050 

6060 

6080 

6100 

6150 

6200 

6200 

6200 

49C40 

28129 

6050 

6060 

6080 

6100 

6110 

6150 

6180 

6190 

48920 

20673 

6033 

6040 

6050 

6060 

6075 

6100 

6130 

6150 

48638 

12587 

6025 

6040 

6050 

6060 

6080 

6100 

6110 

6110 

48575 

4235 

6050 

6050 

6080 

6100 

6100 

6090 

6060 

6060 

43590 

-4237 

6040 

6050 

6075 

6100 

6090 

6050 

6000 

6000 

48405 

-12527 

6040 

6050 

6075 

6075 

6050 

6000 

5980 

5950 

.48220 

-20377 

6030 

6060 

6060 

6050 

6010 

5980 

5900 

5980 

48070 

-27572 

6025 

6030 

6040 

6025 

6010 

5950 

5910 

6000 

47990 

-33933 

6015 

6025 

6000 

6000 

5980 

5930 

5900 

5890 

47740 

-39108 

6005 

6010 

6000 

6000 

5900 

5880 

5880 

5880 

47555 

-43099 

5990 

6000 

6000 

5900 

5890 

5900 

5900 

5900 

47480 

-45860 

5960 

5950 

5930 

5900 

5900 

5940 

5950 

5980 

47510 

-47329 

5950 

5920 

5900 

5900 

5930 

5950 

5990 

6000 

47540 

-47359 

5940 

5920 

5970 

5900 

5950 

5930 

6000 

6080 

47630 

-46005 

5950 

5930 

5910 

5930 

5980 

6000 

6000 

6050 

47750 

-43275 

5950 

5930 

5930 

5950 

5990 

6030 

6040 

6080 

479C0 

-39239 

Table  i;  ;[cor  L,) 

DEFLECTION  STATION  102 
BOUGUER  CORRECTION  COMPON^.NT 


RICE 

.&2£S1_ 

14-21 

14 

15 

16 

17 

18 

19 

20 

21 

SUM 

( 

SUM  X  COS| 

23 

5950 

5950 

5940 

5980 

6000 

6130 

6090 

6190 

48230 

-34103  1 

24 

5960 

5950 

5950 

5990 

6050 

6050 

6150 

6230 

48330 

1 

-2772<:  1 

25 

4 

5960 

5960 

5960 

5990 

6050 

6080 

6090 

6190 

48280 

-20403  I 

1“ 

5980 

5980 

5980 

5990 

6050 

6130 

6180 

6190 

48480 

-12546  j 

5980 

5980 

5990 

5990 

6000 

6100 

6180 

6230 

48450 

-  4224  j 

!  28 

t 

5980 

5990 

6000 

6050 

6010 

6030 

6080 

6:. 00 

48240 

4206  1 

! 

I29 

5980 

6000 

6100 

6100 

6090 

6090 

6050 

6060 

48470 

12544  \ 

|30 

5990 

6000 

6100 

6100 

6100 

6190 

6100 

6160 

48740 

20597  j 

!  3X 

6000 

6000 

6100 

6100 

6150 

6150 

6160 

6230 

48890 

28043  j 

132 

6000 

6000 

6080 

6080 

6150 

6200 

6220 

6260 

48990 

34640  ' 

;33 

6010 

6000 

6030 

6000 

6090 

6190 

6140 

6160 

4862C 

39829  ! 

! 

.  34 

6010 

6010 

6010 

6020 

6030 

6050 

6090 

6140 

48360 

43828  i 

!  35 

6020 

6020 

6030 

6030 

6050 

6075 

6090 

6090 

48405 

46754 

.36 

6030 

6030 

6050 

6050 

6060 

6080 

6100 

6110 

48510 

i 

48325 

I 

I 

I 


1 
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Table  5 

DEFLECTION  STATION  116 
BOUGUER  CORRECTION  COMPONENT 
RICE  14-21 


) 


14 

15 

16 

17 

18 

19. 

20 

21 

SUM  SUM  X  COS  j 

1 

11900 

11450 

11700 

11900 

11500 

10833 

10825 

10425 

90533 

90189 

2 

11900 

11500 

11200 

11200 

11667 

10933 

10425 

10025 

88850 

85820 

3 

12000 

11600 

11250 

11250 

11167 

10767 

10300 

8788 

87122 

78958 

4 

12200 

11700 

11400 

11100 

10917 

10750 

10500 

9650 

88217 

72267 

5 

12500 

12200 

11700 

11150 

10767 

10583 

10550 

9650 

89100 

63002 

6 

12800 

12400 

11900 

11350 

10800 

10367 

10325 

10175 

90117 

51691 

7 

12700 

12550 

12200 

11450 

11067 

10433 

10100 

9975 

90475 

38234 

8 

12450 

12300 

12000 

11600 

11067 

10500 

10162 

9912 

89991 

23289 

9 

12300 

12000 

11700 

11400 

11000 

10533 

10200 

9825 

88958 

7757 

10 

12400 

12050 

11900 

11600 

11167 

10767 

10300 

9838 

90022 

-7850 

11 

12500 

12150 

12000 

11600 

11450 

11000 

10362 

10000 

91062 

-23567 

12 

12100 

11900 

11350 

11538 

11467 

11033 

10600 

10175 

90163 

-38103 

13 

12050 

12300 

11800 

11800 

11767 

10933 

11250 

10825 

92725 

-53187 

14 

12600 

12800 

12400 

12350 

12300 

12167 

11900 

11175 

97692 

-69078 

15 

13300 

13100 

12900 

12600 

12400 

12317 

12012 

11425 

100054 

-81964 

16 

13300 

12850 

12650 

12400 

12250 

11933 

11788 

11575 

98746 

-89493 

17 

13150 

12800 

12600 

12200 

11900 

11816 

11575 

11375 

97416 

-94094 

18 

13100 

12850 

12600 

12300 

12267 

11667 

11362 

11400 

97546 

-97175 

19 

13050 

12700 

12450 

12350 

12100 

11733 

11725 

11900 

98008 

-97635 

20 

13150 

12800 

12500 

12250 

12083 

11917 

11859 

11975 

98534 

-95174 

21 

13150 

12800 

12500 

12350 

12300 

11967 

11744 

11900 

98711 

-89462 

22 

13000 

12700 

12350 

12325 

12133 

11900 

11744 

11662 

97814 

-80129 

-  — ^ 
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Table  5  (cent.) 
DEFLECTION  STATION  116 
BOUCUER  CORRECTION  COMPONENT 


RICE  UnSGS  14-21 

14 

15 

16 

17 

18 

19 

20 

21 

SUM 

SUM  X  COS 

23 

12800 

12600 

12350 

12150 

12083 

11883 

11685 

11950 

97501 

-68943 

24 

12700 

12250 

12175 

12088 

12050 

11900 

11794 

11900 

96857 

-55557 

1 

25 

12500 

12300 

12300 

12175 

12100 

11888 

11850 

11625 

96738 

-40881 

26 

12450 

12500 

12500 

12400 

12300 

11834 

11912 

11925 

97821 

-25316 

27 

12550 

13000 

12800 

12725 

12367 

12000 

11875 

12175 

99492 

-  8675 

28 

12450 

13000 

12900 

12775 

12780 

12200 

12100 

12450 

100655 

8777 

29 

13000 

13000 

12900 

12600 

12500 

12092 

12225 

12512 

100829 

26094 

30 

13050 

13050 

13000 

12600 

12267 

12100 

12325 

12625 

101017 

42690 

31 

13200 

13100 

13360 

12900 

12600 

12467 

12375 

12675 

102677 

.58895 

32 

13145 

13100 

13000 

13150 

12900 

12833 

12375 

12638 

103141 

72931 

33 

13000 

13000 

12700 

12250 

12267 

12500 

12738 

12100 

100555 

82374 

34 

12400 

12500 

12500 

11750 

11367 

11767 

11975 

11825 

96084 

87081 

35 

11850 

11300 

12300 

11750 

11333 

11233 

11250 

11025 

92041 

88902 

36 

11900 

11600 

11550 

13900 

11633 

10950 

11228 

10912 

93673 

93317 

■ „„  ■  » . 
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TABLE  6 
STATION  116 

A0  CORRECTION  -  RINGS  A, 

P 


RING  i  A 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

,  15 

I 

16 

17 

18 
;  19 

20 

21 

i  22 

i 

!  23 


AVG. 
CMPRT . 
HEIGHT 


13250 

13250 

13300 

13450 

13650 

13750 

13850 

13800 

13750 

13700 

13750 

13800 

13800 

13950 

13950 

13950 

13900 

13850 

1  13750 

13750 
13750 
I  13800 
’  13800 


vN  G  B - 

AVG.  AVG. 


AVG.  AVG. 

HEIGHT  CMPRT. 

DIFF.  I  HEIGHT  DIFF. 


.10 

.10 

.15 

.20 

.30 

.30 

.30 

.25 

.25 


12700 
12900 
13100 
13250 
13400 
13550 
13500 
13400 
13150 
13200 
13300 
13400 
13400 
I  13350 
13400 
13550 
13430 
13680 
13600 
13550 
13400 
13450 

I 

i  13400 


DIFF. 

1.35 


B,  AND  C 

^ _ 

AVG. 
CMPRT . 
HEIGHT 

12370 


_<L-C _ 

AVG. 

HEIGHT 

DIFF. 

1.68 


aSp  1 

DEFt.  ; 

CMPMTl2‘ 

.  ...XlO" 

9.66 


1.15  ; 

12430 

1.62  i 

1 

8.21 

.95  1 

12500 

1 

1.55  i 

! 

8.06 

.80 

!  12670 

1.33  1 

3.13 

.65 

I 

i  12930 

1 

1.12 

1.09 

i 

.50 

1 

j  13230 

.82 

.36  i 

1 

.55 

13000 

1.05 

.32  ; 

.65 

12800 

1.25  1 

1 

.34  1 

.90 

12300 

1.25 

.  17 

.85 

12800 

1.25 

.17 

j 

.75 

12970 

1.03  1 

! 

.32  i 

.65 

12570 

1.43  I 

.30  ! 
i  1 

.65 

12430 

1.62  1 

1.33  1 

i 

i 

.70 

12530 

1.52 

j  .  1.36 

.65 

12830 

1.22 

.90 

.50 

13470 

.53 

.20  1 

.62 

1  13500 

.55 

j.  .34 

.37 

i 

13380 

.67 

.26 

.45 

13300 

.75 

.  54 

13300 

13380 

13200 

13100 


i  CJ 
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table  6  (cont.) 
STATION  116 


CORRECTION  -  RINGS  A,  B 

P  _ _ 


‘  TT' 
, 

•  i 

1 

i  i 

R.X-H- 

AVG. 
CMPRT . 
HEIGHT 

G  A  _ 

AVG. 

HEIGHT 

DIFF. 

_ r,.l.b. 

AVG. 

CMPRT. 

HEIGHT 

G  B 

AVG. 

height 

DIFF. 

I _ 1.,. 

1  24  * 

13800 

-25  ■■■"! 

13430 

.62  1 

:a5  i 

13800 

•  25  j 

13430 

.62 

1 

i  26 

13780 

1 

.27  , 

13430 

.62 

1 

'•  27 

13750 

.30 

13400 

.65 

i 

i  28 

13500 

.55 

13380 

.67 

1  29 

13750 

.30 

13450 

.60 

1  30  ! 

13750 

.30  j 

13480 

.57 

j 

1 

13750 

.30 

13480 

.57 

1  32 

j  13800 

.25 

13480 

.57 

> 

1  33 

13700 

.35 

13450 

.60 

i  34 

i  13600 

.45 

1  13100 

.95 

i 

1  35 

13400 

.65 

12800 

1.25 

1 

i  13250 

1 

.80 

12650 

1.40 

,  AND  C. 


^ _ R^I. 

AVG. 
CMPRT . 
height 

N_G _ C _ 

AVG. 

height  1 

DIFF. 

1 

i"  A0p 

1  defl  « 
cmpnt:* 

L— ufl- 

13030 

1.02 

; 

- 

13000 

1.05 

I 

.  .391 

13000 

1.05 

1 

-  •“! 

1 

13000 

1.05 

_  .09  j 

i 

13100 

.95 

I 

.00 

1 

13200 

.85 

.  19 

13280 

.77 

.27 

1  13310 

.74 

.36 

13350 

.70  1 

,  • 

36  1 

13100 

.95 

• 

34 

12870 

1.18 

2, 

.27 

12180 

1.87 

j  7 

.34 

12100 

1.37 

1 

.64 

i 

1 

i 

I 

i 
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table  7 

STATIOiT  116 

CORRECTIOII  -  RINGS  D,  E,  AND  F 


AVG. 

C^^PRT. 

HEIGHT 


11600 

11400 

11800 

12000 

12200 

12000 

12000 

12000 

12000 

12600 

12800 

12800 

12800 

12900 

12700 

12500 

12500 

12700 

13000 

13100 


AVG. 

HEIGHT 

DIFF. 


R  I  N 

G  E 

AVG. 

CMPRT. 

HEIGHT 

AVG. 

HEIGHT 

DIFF. 

11400 

2.65 

10900 

3.15 

10800 

3.25 

.0600 

3.45 

10600 

3.45 

10800 

3.45 

11000 

3.05 

11300 

2.75 

11500 

2.55 

12400 

1.65 

12200 

1.85 

11900 

2.15 

12000 

1 

2,05 

12200 

1 

1.85 

12100 

1.95 

12100 

1.95 

12100 

1.95 

12400 

1.65 

12600 

1.45 

12300 

1.75 

AVG. 

CMPRT. 

HEIGHT 


10400 

10100 

10200 

10100 

10000 

10000 

10200 

10300 

11200 

11800 

11700 

11400 

11800 

11800 

11700 

11700 

11800 

11900 

12200 

12400 


AVG. 

CMPRT. 

DIFF. 


A9p 

DEFL. 

CMPNT 


.109 

.137 

.088 

.066 

.036 

.014 

-.012 

-.034 

-.044 

-.019 

-.019 

-.026 

-.022 

-.017 

-.019 

-.018 

-.012 

-.003 

.002 

.005 


TABLE  7  (cont.) 
STATION  116 


COREECTIOK  -  RINGS  D,  E,  AND  F 


RING  D 

:  R  I  N  G  E 

R  I  N 

G  F 

AVG. 

CMPRT. 

HEIGHT 

AVG. 

HEIGHT 

DIFF. 

AVG. 

CMPRT. 

HEIGHT 

AVG. 

HEIGHT 

DIFF. 

AVG. 

CI4PRT. 

HEIGHT 

AVG, 

HEIGHT 

DIFF. 

AGp 

DEFL. 

CMPNT. 

21 

13100 

.95 

12800 

1.25 

12600 

1.1^5 

1 

.005 

1 

22 

12800 

1.25 

12400 

1.65 

1 

12500 

1.55 

.043 

23 

12200 

1.85 

11300 

2.75 

11500 

2.55 

.054 

1  i 

2U  : 

11900 

2.15  I 

11400 

2.65 

11000 

3.05 

!  .080 

r  1 

III  Mil,  ■ 

m 

O 

o 

o 

O 

o 

in 

o 

o 

o 

1-4 

CO 

CM 

VO 

vO 

cO 

1  CO 

evi 

in 

o 

in 

o 

1-4 

eo 

o 

St 

CO 

CM 

in 

f-4 

00 

f  H  ■ 

»o 

m 

m 

m 

o 

VO 

m 

xd- 

in 

ex 

Cv 

o 

CO 

f-4 

VO 

VO 

1  • 

v-4 

1-4 

t-4 

1-4 

1 

f-4 

1  I 

o 

o 

o 

o 

o 

o 

o 

in 

m 

m 

r4 

CO 

CO 

CO 

CO 

CO 

1  rsi 

»o 

m 

o 

o 

o 

o 

r«. 

1-4 

rx 

rx 

CO 

in 

CO 

rx 

CO 

1-4 

H  ! 

m 

m 

in 

m 

in 

in 

m 

m 

St 

VO 

CO 

ON 

r-4 

St 

VO 

St 

1 

f-4 

f-4 

f-4 

i 

!  1 

o 

o 

o 

o 

o 

o 

in 

m 

•n 

m 

CO 

CO 

CM 

o 

rx 

m 

[VO  . 

m 

m 

m 

in 

o 

CM 

CM 

CO 

o 

CO 

CO 

CM 

o 

f-4 

St 

H  1 

»o 

m 

m 

m 

m 

in 

in 

in 

St 

m 

CO 

o 

CM 

in 

OV 

CO 

1 

1-4 

1-4 

1-4 

1-4 

!  ! 

1-4 

1 

1 

o 

m 

m 

o 

o 

o 

m 

m 

o 

o 

VO 

•-4 

o 

rx 

rx 

1-4 

:  m  i 

lO 

CM 

CM 

o 

o 

o 

o 
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116 

117 

117 

113 

23 

111 

112 

112 

113 

113 

113 

24 

106 

107 

107 

107 

10  6 

107 

25 

103 

102 

102 

101 

100 

99 

26 

100 

97 

96 

93 

92 

90 

27 

55 

52 

90 

39 

0  n 

J  J 

36 

23 

90 

39 

3  6 

34 

32 

76 

29 

n  n 

36 

33 

81 

75 

70 

30 

86 

34 

31 

76 

63 

65 

31 

35 

33 

77 

71 

65 

59 

32 

34 

31 

76 

69 

63 

56 

33 

34 

31 

76 

70 

62 

5  6 

34 

3  4 

31 

77 

73 

66 

53 

35 

r>  r 

U  J 

82 

30 

75 

71 

60 

36 

0  D 

33 

32 

79 

75 

65 
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3  5 

Sum  X  Sln4 

64 

111.44 

75 

337.43 

31 

560.79 

34 

77U.92 

33 

979.33 

92 

1163.26 

1329.54 

110 

1476.36 

123 

1589.94 

129 

1634.76 

132 

1612.09 

135 

1537.99 

139 

1409.84 

146 

1237.43 

153 

1021.01 

160 

762.79 

160  469.20  i 


153  157.05  I 

146  -  154.95  i 

136  -  452.12  ! 

129  -  724.76  i 

123  -  965.37  | 

115  .1160.35  ! 


107  .1300.07  , 

100  .1383.45  ■ 


^  n 

^  sj 

.1419.87 

31 

.1416.60 

72 

..  1353.81 

63 

.  1275.95 

60 

.1156.44 

53 

. 1014.99 

50 

.  861.25 

50 

.  695.78 

52 

.  516.84 

55 

-  320.65 

56 

.  109.78 

page  25 


Table  13 

DEFLECTIOl!  STATIOI^ 

DEFLECTIOIl  r.ESULTS  niCE  RIMGS  22-29 


Complete  Bou3uer  Anomaly  Component 


2  2 

23 

24 

25 

26 

27 

23 

29 

1 

120 

113 

116 

115 

113 

110 

105 

102 

2 

120 

120 

113 

116 

115 

112 

111 

106 

3 

121 

121 

120 

113 

118 

114 

113 

111 

k 

122 

122 

122 

121 

120 

117 

115 

113 

5 

123 

123 

123 

123 

122 

121 

119 

116 

6 

125 

125 

125 

125 

125 

125 

123 

121 

7 

126 

120 

127 

127 

120 

120 

123 

127 

8 

127 

123 

129 

130 

131 

131 

131 

131 

9 

127 

129 

130 

131 

131 

133 

133 

134 

10 

123 

129 

131 

131 

132 

133 

135 

136 

11 

129 

130 

131 

131 

133 

134 

136 

133 

12 

129 

150 

131 

131 

133 

135 

137 

140 

13 

129 

129 

131 

131 

133 

135 

130 

140 

14 

129 

129 

130 

131 

132 

135 

130 

140 

15 

129 

129 

129 

130 

132 

134 

137 

140 

16 

123 

129 

129 

130 

131 

133 

135 

130 

17 

127 

120 

129 

129 

129 

132 

134 

136 

13 

127 

127 

123 

129 

129 

130 

132 

133 

19 

127 

127 

127 

120 

120 

129 

129 

131 

20 

126 

126 

12  6 

127 

127 

123 

128 

129 

21 

125 

125 

125 

125 

125 

12  6 

127 

127 

22 

124 

124 

124 

124 

124 

124 

124 

125 

23 

123 

123 

123 

123 

123 

123 

123 

122 

24 

122 

122 

122 

122 

121 

121 

121 

120 

25 

122 

121 

121 

120 

120 

119 

119 

110 

26 

121 

121 

120 

119 

113 

117 

116 

115 

27 

121 

120 

119 

118 

116 

115 

115 

114 

28 

120 

119 

117 

116 

115 

114 

113 

112 

29 

119 

113 

116 

115 

114 

113 

112 

110 

30 

113 

116 

115 

114 

113 

111 

110 

107 

31 

117 

116 

115 

113 

111 

109 

107 

105 

32 

117 

115 

114 

112 

110 

103 

106 

103 

33 

117 

115 

114 

111 

109 

107 

104 

100 

34 

117 

116 

114 

111 

109 

106 

102 

99 

35 

113 

116 

115 

112 

110 

107 

101 

99 

36 

119 

117 

115 

113 

111 

103 

103 

99 

% 


Table  13  (cpnt.) 
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DEFLECTION 

Complete 

RESULTS  RICE  RINGS 

Bouguer  Anomaly  Compo 

33  34 

30-3  7 

nent 

35 

36 

30 

31 

32 

Sum  X  61n( 

1 

101 

100 

98 

94 

90 

37 

33 

-  135.33 

2 

103 

102 

98 

95 

92 

r»  n 

U  ' J 

35 

409.^16 

3 

105 

103 

99 

97 

93 

39 

36 

679.54 

4 

112 

106 

102 

98 

95 

92 

88 

943.57 

5 

114 

112 

106 

102 

93 

95 

92 

1194,29 

6 

113 

115 

112 

107 

104 

100 

97 

1431.14 

7 

125 

122 

117 

114 

111 

107 

104 

1730.28 

3 

131 

130 

127 

123 

119 

116 

113 

1832.31 

9 

135 

135 

135 

133 

130 

126 

124 

1958.53 

10 

137 

139 

140 

140 

138 

133 

136 

2015.31 

11 

140 

141 

143 

147 

143 

150 

152 

2011^97  ! 

12 

141 

144 

145 

153 

155 

159 

152 

1925.89 

13 

142 

145 

147 

157 

161 

165 

170 

1763.74 

14 

142 

145 

147 

159 

165 

170 

175 

1532.29 

15 

142 

144 

146 

159 

167 

173 

179 

1244. 71 

16 

141 

143 

145 

157 

165 

173 

100 

911.55 

17 

140 

142 

144 

155 

162 

171 

131 

553.57 

18 

136 

139 

142 

150 

157 

170 

131 

103.99 

19 

133 

136 

140 

145 

153 

165 

173 

-  176,67 

20 

130 

133 

136 

139 

146 

157 

171 

-  525.11 

21 

123 

129 

131 

136 

140 

149 

162 

-  836.75 

22 

126 

126 

128 

129 

134 

142 

147 

-1104.18 

23 

122 

122 

123 

123 

125 

132 

140 

-  1322.98 

24 

120 

119 

113 

119 

120 

124 

127 

-  1439^31 

25 

117 

116 

116 

116 

116 

117 

120 

-  1611.40 

26 

115 

114 

113 

113 

113 

113 

113 

-1681.63 

27 

113 

112 

110 

107 

106 

105 

105 

-  1639.  56 

28 

111 

110 

106 

103 

99 

98 

96 

-  1642.  73 

29 

109 

106 

91 

97 

94 

92 

88 

-  1539.65 

30 

105 

103 

33 

94 

90 

33 

35 

-1411.11 

31 

103 

100 

36 

93 

39 

3  5 

33 

-  1255.01 

32 

100 

97 

34 

91 

36 

81 

77 

-  1061.36 

33 

90 

95 

90 

38 

35 

84 

79 

-  850.11  ! 

34 

97 

94 

91 

33 

3  7 

85 

82 

-  633.06  j 

35 

97 

95 

93 

90 

88 

35 

83 

-  390.53  ; 

36 

98 

97 

95 

92 

89 

36 

32 

-  132.89 
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DEFLECTION  STATION  2 
BOUGUER  CORRECTION  COMPOHEiTT 
RICE  MNGS  15-21 


15 

16 

17 

18 

19 

20 

21 

SUM 

SUM  X 

1 

720 

720 

770 

800 

748 

730 

695 

5133 

452 

2 

730 

710 

750 

750 

723 

725 

690 

5033 

1315 

3 

750 

700 

720 

750 

770 

750 

730 

5170 

2135 

4 

7G0 

720 

700 

740 

776 

770 

755 

5221 

2995 

5 

760 

730 

720 

720 

770 

742 

700 

5142 

3636 

6 

740 

730 

700 

690 

715 

725 

662 

5042 

4130 

7 

740 

700 

700 

660 

652 

672 

640 

4764 

4313 

3 

720 

710 

750 

700 

670 

740 

633 

432J 

4663 

5 

7  00 

750 

300 

750 

725 

650 

625 

5060 

5041 

10 

770 

750 

750 

650 

660 

640 

505 

4343 

4325 

11 

730 

710 

650 

600 

605 

642 

610 

4627 

4469 

12 

570 

650 

650 

660 

700 

535 

655 

4630 

4241 

13 

670 

650 

670 

750 

770 

715 

690 

4915 

4026 

14 

670 

630 

700 

700 

750 

732 

742 

4974 

3517 

15 

700 

740 

720 

730 

730 

7  65 

775 

5160 

2960 

16 

720 

770 

700 

300 

GOO 

730 

303 

5463 

2311 

17 

640 

770 

300 

050 

860 

r»  ri  o 
£, 

900 

5702 

1476 

10 

650 

750 

750 

300 

365 

950 

950 

5715 

493 

19 

630 

750 

750 

760 

73  2 

332 

812 

5365 

-  463 

20 

700 

700 

750 

750 

762 

750 

732 

5194 

-1344 

21 

730 

700 

740 

300 

330 

920 

900 

5670 

-2325 

22 

760 

750 

730 

700 

730 

340 

903 

5543 

-3132 

23 

730 

24 

730 

25 

730 

2C 

770 

27 

790 

23 

820 

2C 

330 

30 

050 

31 

350 

32 

340 

33 

310 

34 

780 

35 

740 

36 

730 
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Table  Ih  (cont.) 
di:fl2Ctioi:  statioh  2 
BOUIUER  CORR’UCTIOi!  COMPOUEiTT 


RICE 

'RINOS  15- 

21 

16 

17 

13 

ID 

20 

21 

sun 

SUM  X  Sill 

310 

760 

760 

790 

343 

925 

5673 

-4011 

320 

340 

7  60 

755 

730 

3  63 

5  603 

-4590 

760 

3  60 

330 

755 

792 

912 

5639 

-5111 

750 

750 

310 

330 

770 

325 

5505 

-5317 

300 

300 

300 

790 

735 

040 

5505 

-5534 

340 

850 

350 

370 

742 

0  63 

5040 

-5313 

350 

3  50 

350 

530 

345 

335 

5390 

-  5639 

860 

300 

340 

315 

352 

303 

5325 

-  5279 

790 

7  60 

3  50 

320 

7  50 

740 

5570 

-4563 

760 

750 

800 

858 

775 

730 

5513 

-3393 

770 

740 

010 

3  10 

760 

735 

5435 

-3146 

730 

750 

840 

010 

760 

753 

5423 

-2294 

730 

750 

310 

3  03 

780 

760 

5373 

-1392 

725 

7  60 

830 

065 

750 

745 

5405 

-  471 

